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Abstract. The two-point function for tensor metric perturbations around de Sitter spacetime 
including one-loop corrections from massless conformally coupled scalar fields is calculated 
exactly. We work in the Poincare patch (with spatially flat sections) and employ dimensional 
regularization for the renormalization process. Unlike previous studies we obtain the result 
for arbitrary time separations rather than just equal times. Moreover, in contrast to existing 
results for tensor perturbations, ours is manifestly invariant with respect to the subgroup of 
de Sitter isometrics corresponding to a simultaneous time translation and rescaling of the 
spatial coordinates. Having selected the right initial state for the interacting theory via an 
appropriate ie prescription is crucial for that. Finally, we show that although the two-point 
function is a well-defined spacetime distribution, the equal-time limit of its spatial Fourier 
transform is divergent. Therefore, contrary to the well-defined distribution for arbitrary 
time separations, the power spectrum is strictly speaking ill-defined when loop corrections 
are included. 
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1 Introduction 

The amplification of quantum vacuum fluctuations in inflationary models provides a natural 
mechanism for the generation of primordial cosmological inhomogeneities. This constitutes 
the key element of the currently most successful paradigm for explaining the observed cosmic 
microwave background (CMB) anisotropics and large scale structure of the universe. The 
standard theoretical analysis and calculation of the primordial spectrum relies on the rela- 
tively simple procedure of quantizing linear cosmological perturbations, which amounts to a 
tree-level calculation in an accelerating FLRW background (typically with quasi-exponential 
expansion and, hence, close to de Sitter spacetime). This framework also predicts quite 
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generically the existence of an approximately scale-invariant background of primordial gravi- 
tational waves (tensor perturbations). The situation, however, becomes much more complex 
and subtle when quantum effects at nonlinear level are taken into account and loop cor- 
rections are considered. In fact, a clear understanding of how significant the effects in this 
context are is still an unsettled issue. 

Tsamis and Woodard proposed that radiative corrections from graviton loops at two- 
loop order and higher could generate a secular screening of the cosmological constant [1, 2]. 
The basic idea is that in a certain gauge free gravitons behave like two massless minimally 
coupled scalar fields (one for each polarization), and more and more modes get excited as 
they leave the horizon without being diluted away by the exponential expansion. It is then 
argued that due to the nonlinear nature of the gravitational interaction, the back-reaction 
of the attractive gravitational field between the growing number of IR graviton modes slows 
down the accelerated expansion [3]. 

In addition, there has recently been renewed interest in the possibility that quantum 
effects in interacting theories could lead to an instability even for massive theories evolving on 
a fixed de Sitter background [4]. This possibility has been analyzed both at tree level [5-10] 
and including loop corrections [11, 12]. In this context, the use of Euclidean methods has 
proven very useful. They provide a natural and direct way of generating the Hartle-Hawking 
state [13] for the interacting theory, which generalizes the de Sitter-invariant Bunch-Davies 
vacuum of the free theory to the interacting case. It exhibits de Sitter-invariant quantum 
correlation functions, which can be obtained by analytic continuation from Euclidean to 
Lorentzian time, and has been shown to constitute a late-time attractor (defined in a precise 
sense) for the evolution of generic initial states [14-16]. Although these studies focused on 
sufficiently massive theories (depending on the coupling strength) , generalizations to massless 
and very light fields have also been developed [17, 18]. 

Our goal is to extend these considerations to the gravitational case. As a first step, we 
will concentrate on the effects of matter loops (neglecting graviton loops) on the two-point 
function characterizing the quantum fluctuations of the metric perturbations around a de 
Sitter background. There are clearly significant differences between the main case considered 
in the studies mentioned in the previous paragraph, a scalar field with no derivative interac- 
tion, and the gravitational case. One is the (partial) derivative character of the gravitational 
interaction, which improves for instance the IR behavior of massless minimally coupled mat- 
ter fields in this case. (In addition, perturbative quantum gravity around a given background 
is power-counting non-renormalizable and should be treated as a low energy effective field 
theory [19, 20], but this aspect is shared by those analyses that considered scalar fields with 
nonlinear interactions in arbitrary dimensions [21].) A second difference is the existence of a 
gauge symmetry associated with diffeomorphism invariance and the need to consider appro- 
priate gauge-invariant observables, which is a rather nontrivial aspect even in perturbative 
quantum gravity [22]. Moreover, one should restrict one's attention to "sufficiently local" 
observables that properly characterize the geometrical properties within a region of finite 
physical size. This point has been crucial in order to construct IR-safe observables in situa- 
tions which would otherwise lead to divergences in the absence of an IR cut-off [22-25]. Here 
we calculate the one- loop correction to the tensorial metric perturbations, which is the key 
ingredient to obtain such kind of observables. The scalar and vectorial metric perturbations, 
in contrast, can be directly obtained from the stress tensor correlation function. 

Our strategy (mainly for computational simplicity) will be to do the calculation using 
spatially flat coordinates in the Poincare patch of de Sitter spacetime. We will employ the 
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closed-time-path (CTP, also known as in-in) formalism, which gives the real time evolution of 
true expectation values and correlation functions rather than transition matrix elements [26] . 
Moreover, an asymptotic initial state selected by an appropriate Wick-rotation prescription 
for the two CTP branches will be considered. This corresponds to an asymptotic adiabatic 
vacuum of the interacting theory which is expected to coincide with the Hartle-Hawking state 
generated by a Euclidean path integral. In fact, the equivalence between such a Euclidean 
calculation of the correlation functions (followed by an analytic continuation to Lorentzian 
time) and the in-in calculation has been established for massive scalar theories in [21]. That 
result, however, cannot be directly extended to the gravitational case because the graviton 
propagator does not fall-off quickly enough for large spacetime separations. Nevertheless, 
we expect that as long as one considers gauge-invariant observables sensitive only to geo- 
metric properties localized within a region of finite physical size, the two methods will yield 
equivalent results. 

The effects of loops of conformal fields on tensor cosmological perturbations have re- 
cently been studied in [27], where a result was found which grows as the initial time tends 
to minus infinity and diverges in that limit. That kind of result together with the small am- 
plitude of the observed CMB anisotropics would impose a bound on the possible number of 
e-foldings during inflation. It is also incompatible with de Sitter invariance, which implies (as 
a necessary but not sufficient condition) that there should be a (co-moving) time translation 
symmetry for spatially flat sections provided that one compares properties at different times 
involving the same physical distance and momentum scales. 

Exact de Sitter invariance will only happen for a very particular choice of initial state 
corresponding to the de Sitter-invariant vacuum of the interacting theory (rather than the 
free one): in general results will not be exactly invariant for any other state. However, if such 
a de Sitter-invariant state exists and is a late-time attractor, in analogy with the results for 
scalar theories [15, 16] and as suggested by the intuition that regular initial excitations get 
red-shifted away by the exponential expansion, one expects that for sufficiently regular initial 
states the result for observables characterizing the local geometry in a finite region should 
become insensitive to the initial conditions and tend at sufficiently late times (or when the 
initial time tends to minus infinity) to the de Sitter-invariant result. 

These are of course just expectations which should be confirmed by an actual compu- 
tation. Here we provide an explicit and exact calculation (at one loop in the matter fields) 
with a result which is well behaved in the limit of asymptotic initial time and respects the 
translation symmetry (rescaling symmetry for co-moving momenta and conformal times) 
mentioned above. Our result also confirms the absence of terms involving a logarithmic 
dependence with respect to the co-moving momentum only. (Such terms, which violate the 
rescaling symmetry, were originally found in calculations of one-loop corrections to the power 
spectrum of scalar [28, 29] and tensor [30] cosmological perturbations, but were later argued 
not to be correct in [31], where their absence was explicitly shown for scalar perturbations.) 
Furthermore, we will explain in detail how the behavior found in [27] is a consequence of not 
considering the full contribution of matter fields to the quantum correlation function of the 
metric perturbations (at one loop). 

It is worth emphasizing that we have obtained the correlation function for two arbitrary 
times rather than just the equal-time limit, as done in previous studies. This is an important 
point if one wants to check full de Sitter invariance: by considering only equal times, instead, 
one can only check invariance under a subgroup of the de Sitter group. Moreover, together 
with an explicit implementation of the renormalization procedure, our exact result for two 
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arbitrary times has enabled us to uncover a previously overlooked subtlety: although the 
one-loop two-point function is a well-defined spacetime distribution, the power spectrum 
(the equal-time limit of its spatial Fourier transform) is divergent, even though all usual 
divergences arising in the calculation have been properly regulated and canceled out by local 
counterterms. 

The paper is organized as follows. In section 2 we briefly review the in-in (or CTP) 
formalism and the ie prescription which selects in this context the adiabatic vacuum of the 
interacting theory as the asymptotic initial state. Next, a general expression for the quantum 
two-point function at one-loop is derived in this framework. It can be conveniently written 
in terms of the effective action that results from integrating out the matter fields. Such an 
effective action is presented in section 3 for metric perturbations around a FLRW background 
when integrating out a conformal scalar field. Making use of these tools, in section 4 we 
compute the two-point function of the tensor perturbations around a de Sitter background 
for an arbitrary pair of points and including the one-loop correction from a conformal field. 
Moreover, we give a detailed explanation for the discrepancy of our results as compared to 
those of [27]. Finally, we discuss our main findings in section 5. A number of technical 
aspects are contained in the appendices. Furthermore, a thorough discussion of the fact 
that the power spectrum is, strictly speaking, ill defined when including loop corrections is 
provided in appendix F. 

We use the "-|--|— 1-" sign convention of [32], but use only Latin tensor indices which 
range over space and time. Throughout the paper we work in natural units c = h = 1 and 
take = IGttGn, where Gn is Newton's constant. 



2 A general overview of the formalism and calculation 
2.1 The "in-in" formalism and the ie prescription 

In the quantum field theoretical treatment of scattering problems one usually calculates the 
transition matrix element of an operator between two states, 

{aim, (2.1) 

where \a) and can be taken to be two different in and out vacuum states, with the 
particle content of the real states incorporated into the operator A. Up to a (possibly 
infinite) proportionality constant, one has the path integral representation 

(out|i|in) ~ yA[0]e'^['^]p0, (2.2) 

assuming A to be time-ordered (which is no problem in practice) and denoting by (p the field 
content of the theory, and by S the action. 

If y4 is a polynomial in the fields (and their derivatives), or can be well approximated by 
a polynomial, one can calculate its transition matrix elements by adding a classical source J 
to the action 

S^S + j J{x)4>{x)d'^x, (2.3) 
where n is the number of spacetime dimensions, and constructing the generating functional 

Z[J] = [ e'^['^'-^]p</., (2.4) 
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in terms of the external source . Normalized matrix elements of time-ordered polynomials 
in the fields are then obtained by functionally differentiating with respect to the source J, 



(out|T(A(xi) • • • 4>{xn)\in) 



(2.5) 



j=o 



(out I in) 

so that we can also display the generating functional as 

Z[J] = (out|re^^'^(^)'^(^')'^^|in) . (2.6) 

To implement the standard flat-space choice of vacuum in the path integral (2.2) or 
(2.4), one slightly tilts the time integration contour on the complex plane to include an 
imaginary part 

t^t(l-ie) (2.7) 

with e > (see for instance section 4.2 in [33]). This selects the asymptotic vacuum as the 
state of lowest energy of the full interacting theory, which includes appropriate correlations 
between the different fields or even different modes of the same field. ^ In practice, one 
calculates the integral from some initial time to to some final time T, and takes the limits 
to — — oo, T — )• oo in a slightly imaginary direction; or, as it is often done, rotate the time 
axis altogether onto the imaginary axis, which gives rise to Euclidean quantum field theory. 

In a cosmological setting, one is interested in true correlation functions, which are ex- 
pectation values of operators rather than transition matrix elements. Moreover, one typically 
needs to impose initial conditions at early times instead of boundary conditions at both early 
and late times. Furthermore, as was explicitly shown in [4, 7] for de Sitter space, in an 
exponentially expanding spacetime "in-out" perturbation theory has an infra-red divergence 
because of the expanding space volume. For all these reasons we are naturally led to consider 
the "in-in" formalism, where one specifies initial conditions at some initial time (in certain 
cases, such as the exponentially expanding patch of de Sitter with spatially flat sections, one 
can specify these initial conditions at past infinity and define an asymptotic "in" vacuum). 

In order to calculate an expectation value using path integrals, one inserts the identity 
operator as a sum over an orthonormal basis of states at some "final" time T, 



(in|j4|in) = (in|a, T){a, r|j4|in) = ((a, T|in))* {a, T\A\m) 

a a 

~ j ^[(/>+](5[(/>+(r) -a](5[0^(T) -a]e^^['^^le-'^[<^"]p0+P(/)-Pa, 



(2. 



where the states \a,T) correspond to an orthonormal basis of field-configuration eigenvectors 
in the Heisenberg picture, such that (l){T)\a,T) = a\a,T). Since we have two path integrals 
for each degree of freedom, we need two copies of the fields which we have labeled (j)'^ and 
<j)~ . Instead of enforcing the separate equality of both fields to a and then integrating over 
all field configurations a at time T, in the following we can directly enforce the equality of 
the fields d)~^ and (h~ at that time. 



'^In a gravitational context, general covariance of the action (2.3) requires that the source J{x) transforms 
as a density under general coordinate changes. 

^This procedure selects the ground state of the interacting theory only for a time-independent Hamiltonian. 
Under the appropriate conditions, however, in the time-dependent case, this prescription can still select an 
adiabatic vacuum of the theory at early times. 
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Similarly to the "in-out" case, a generating functional for the correlation functions can 
be introduced, but one needs two sources and J~ corresponding to each one of the two 
path integrals. In addition, because of the complex conjugation of the second path integral's 
integrand, which effectively reverses the time integration in the action, the corresponding 
operator expression is naturally anti-time-ordered, which we denote by the symbol T. Fur- 
thermore we may generalize the expression by considering a density matrix p as the initial 
state, which for an initially pure state reduces to |in)(in|. Altogether, the corresponding 
"in-in" generating functional is given by 




= j 4,/.+ (r)-r(T)]4</.+ (to)-0^]4'^"(*o)-0o]W|p|0o>x (2.9) 

where 4>^ are the field configurations at the initial time to- This time can be finite, although 
one would need in principle to consider an appropriate dressed state for the interacting theory. 
In this paper instead we will take to = — oo and use an appropriate prescription to select the 
right asymptotic initial state, just as in the case described above. 

The expression in equation (2.9) can be interpreted as follows: starting from an initial 
state described by the density matrix p, we evolve it forward in time under the influence of 
an external source J+ up to some flnal time T, and then evolve backwards in time under 
the influence of a different source J~ (were the two sources the same, we would have clearly 
achieved nothing). Because of this, the "in-in" formalism is also known as the closed-time- 
path (CTP) formalism, where the usual time integration contour going from — oo to +00 is 
replaced by one going from —00 to T and turning back to —00. In fact, one can alternatively 
write the path integral expression for the CTP generating functional in equation (2.9) in 
terms of a single copy of the fields by having the time integrals in the action run over this 
whole CTP contour (and letting the fields and the currents have independent values on both 
halves of the contour) [34] . Similarly, since Feynman rules can be directly derived from the 
generating functional, one can show that the Feynman-Stiickelberg diagrams for the "in-in" 
formalism can be simply obtained from the "in-out" diagrams by evaluating the time integrals 
associated with the interaction vertices along the full CTP contour. 

Functionally differentiating the CTP generating functional Z[J^ ^ J^] with respect to 
iJ^ and — iJ~ one now generates expectation values of ordered polynomials in the fields, 
which in general are path- ordered (denoted by V) along the CTP contour rather than time- 
ordered. For the 0"^ fields, the path ordering is the standard time ordering, while for the 
fields it is anti-time ordering as discussed above; and </>~ fields are always to be ordered as if 
they occurred at a later time than any (j)'^ field (see figure 1). 

The ie prescription can also be carried over to the CTP formalism. However, due to the 
complex conjugation of the integrand in the path integral for the 0~ fields, we need to take 
the complex conjugate prescription for them 

t^t(l + ie). (2.10) 

One therefore has to integrate along a contour going from t^ to T, returning to tp with 
tg = (t^)*, and taking at the end of the calculation — — oo(l — ie). The dependence on 
T disappears in the final result as long as it is larger than all times of interest (i.e., all the 
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times in the arguments of the correlation fmiction one wants to calculate), as required by 
causality. The deformed contour and the path ordering are shown in figure 1. 

Note that this prescription is suitable for selecting the asymptotic initial state in space- 
times where the behavior of the modes for free fields is dominated at past infinity by the same 
kind of oscillatory behavior as in Minkowski space (factors with a power-law or weaker time 
dependence are allowed); it then selects the adiabatic vacuum for the interacting theory as 
t — >• — oo. This is the case for the exponentially expanding patch of de Sitter spacetime with 
spatially flat sections, but would not be appropriate for global de Sitter, which is instead 
exponentially contracting coming from past infinity. 

For further details about the CTP formalism see for instance [26, 34-38]. 

^0 0+(r2) 



3 



0"(r4) (r-i) 



Figure 1. The deformed CTP integration contour, together with 4 path-ordered fields. For any 
permutation of the indices we have V(j)'^{Ta)4>'^{Tb)(f)'^{Tc)4>'^{Td) = {T4)(j)~ {T3)(j)'^ {t2)4>'^ (ti) . Here, 
(/)"'" (tq) and ^~(T(j) simply represent the operator (j) evaluated at a time lying, respectively, on the 
upper and lower branch of the deformed complex contour for time integration, and correspond to the 
fields and respectively in the path integrals. 



2.2 Calculating the two-point function at one-loop order 

Taking a mean field approach to quantum gravity, we will quantize the metric perturbations 
around a semiclassical background including the effects of matter loops. Although perturba- 
tive quantum gravity is power-counting nonrenormalizable, it can be consistently studied as 
an effective field theory (EFT) describing quantum gravitational phenomena with character- 
istic length scales well above the Planck length [19, 20]. In order to do so, one needs to 
introduce local counterterms in the bare action for the metric and the matter fields with an 
arbitrary number of derivatives and powers of the curvature, suppressed by the corresponding 
negative power of the Planck mass (or positive power of the Planck length) . The key point is 
that for phenomena with a characteristic length scale L only a finite number of counterterms 
needs to be considered to achieve a certain precision, roughly given by (/p/L)^" with n being 
the total number of pairs of derivatives and/or powers of the curvature. 

More specifically, we want to calculate the two-point function of the metric perturbations 
around a fixed background g'^^^ (so that the total metric is gab = gab + /^^ab) coupled to a 
scalar field. Therefore, we will consider the following CTP generating functional containing 
external sources for the metric perturbations hab (since we are not interested in the matter 
field correlations, we do not need sources for the matter field (J)): 

Z[J+,J-] = 1 6[<p-^iT)-riT)]5[h+{T)-h-iT)]5[h^{to)-h^]6[<P^ito)-4]x 

(2.11) 
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where S[h, (j)] denotes the total action for metric perturbations and matter fields, 4>^ and 
correspond to the field configurations at the initial time to, and p is the density matrix 
specifying the initial state. Note that for simplicity we adopt throughout this subsection 
a shorthand notation with tensor indices, spatial dependence of the fields and spacetime 
integrals omitted. Our calculation will be performed in two steps: first, we explicitly integrate 
over the matter fields (/>^, and next we perform the remaining path integrals for the metric 
perturbations h^. It is easier to proceed if one assumes that the initial density matrix 
factorizes so that ( 0,| , | p | (/>q , /iq ) = ( '/'o^ | P,/, | ) ( ^o" | P/i | ) • We are actually interested in 
states of the interacting theory which involve correlations (entanglement) between the matter 
fields and the metric perturbations, for which the density matrix does not factorize. However, 
if one employs the method described in the previous subsection to select the adiabatic vacuum 
of the interacting theory, one can consider a factorized initial state involving the product of 
the Bunch-Davies vacua for the matter fields and the metric perturbations both treated as 
free fields, integrate out the matter fields more easily and only at the final stage deform the 
CTP integration contour (as shown in figure 1) for the vertices appearing in the Feynman- 
Stiickelberg diagrams associated with the perturbative calculation of the path integral of the 
metric perturbations. In doing so we will neglect contributions from graviton loops. This 
can be formally implemented in a natural way by considering a large expansion for 
matter fields interacting with the gravitational field [39, 40]; for instance, the lowest-order 
contributions to the connected two-point function of the metric perturbations are of order 
whereas any contribution including graviton loops is suppressed by higher powers of 

1/N. 

After integrating out the matter fields, which amounts to a one-loop calculation, one 
is left with the following expression for the CTP generating functional (dropping an infi- 
nite proportionality constant, which anyway will cancel in the calculation of the correlation 
function): 



Z[J+, J-] = J 6[h+{T) - h-{T)]5[h^{to) - 



(2.12) 



where 5*0 [^] is the gravitational part of the original bare action S[h,<j)], i.e. all the terms 
depending only on the metric, and E is the result of the functional integration over the 
matter fields, truncated at second order in h. In the course of "integrating out" the matter 
fields (j), the ultraviolet (UV) divergences associated with matter loops which appear in U 
must be regularized and appropriate counterterms must be included in the bare gravitational 
action Sg- To leading order in 1/A^, it is only necessary to introduce counterterms in 6*0 [/i] 
at most quadratic in the curvature^, so that the bare gravitational action is schematically 
given by 



ao 



j{R- 2/1)7=^ d"x + a2K^ J i?2^/=5d"x + O {r^ 



(2.13) 



where at are dimensionless bare parameters and has been used as a shorthand for terms 
quadratic in the Riemann tensor with all possible contractions. We have introduced a rescaled 



^Since perturbative quantum gravity around a given background is power-counting non-renormalizable, for 
each additional term in the expansion one would need to include counterterms with higher and higher 
powers of the curvature and number of covariant derivatives in Sci^]- Moreover, in that case one would also 
need to include in the full bare action S[h,(l)] counterterms coupling h and <j!>, and with a higher and higher 
number of covariant derivatives. 
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gravitational coupling constant 

R^ = Nk^, (2.14) 

which helps to organize conveniently the calculation when using a expansion.'^ Having 
done so, the divergent terms cancel out in the sum 

5eff[/i+, h-\ = SG[h+] - SG[h-] + S[h+, h~\ , (2.15) 

and the regulator can be taken to its physical value. Note that X! provides an effective 
propagator for h which includes the interaction with 0, so that S^s can be called a CTP 
effective action. 

From now on we will only consider Se^lh'^ ,h^] through quadratic order in the metric 
perturbations since higher-order terms do not contribute to the generating functional and 
the connected correlation functions to leading order in 1/A^, as mentioned above. To analyze 
more closely the different contributions to S'eff, let us expand it in powers of /i, 



(2.16) 



where the superscript denotes how many powers of h appear. The sum of all the first 
order terms corresponds to integrating and h~ with the functional derivative of the CTP 
effective action for the matter fields in a semiclassical background with the metric of the two 
branches equal to (7^^ after differentiation. When equated to zero, such a functional derivative 
corresponds to the semiclassical Einstein equation which governs the quantum back-reaction 
of the matter fields on the mean background geometry [44, 45] 



^2 



Ga,m-Agl, = ^lf^,[g']) , (2.17) 



where the stress tensor expectation value comes from the functional derivative of E with 
respect to hab and here also includes any finite contributions from the counterterms in 
other than the Einstein tensor or the cosmological constant. We will always consider a 
background g^^^ which is a solution of the semiclassical equation (2.17). Therefore, the sum 
of the first-order terms on the right-hand side of equation (2.16) will vanish. Such a choice 
will also guarantee that after the functional integration there will be no term linear in J in 
the exponent of the CTP generating functional, which will be purely quadratic at the order 
at which we are working. 

Furthermore, the zeroth-order terms in equation (2.16) can be factored out of the path 
integral (2.12) and give a factor independent of J which does not contribute to the corre- 
lation functions, obtained by functionally differentiating with respect to J. Disregarding 
them we are, thus, left with a purely quadratic expression for S'cff. Given an exponent 
quadratic in the perturbations and h~ , one can easily perform the functional integration 
in equation (2.12). In order to do so, it is convenient to adopt a matrix formulation with 

hA = {h'^,h^), Ja = {J^,—J~) (note the minus sign), Amn = ( ^ i ^ j and sum over 

*Note that having a fundamental length-scale R characterizing the breakdown of the low-energy gravita- 
tional EFT which is much larger than the effective Planck length Ip = y/T&TTK determined at low energies can 
happen naturally in braneworld models with warped extra dimensions [41] or models with a large number of 
fields [42]. In fact, in some cases both possibilities can be understood as equivalent descriptions within the 
framework of the AdS/CFT correspondence [43]. 
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repeated capital indices regardless of their position. Expanding now the effective action S^s 
in powers of k, 

Ses[h~^,h-] = - S^Q\h-] + U^^^^h+,h'] = So[h+,h-] + K^S2[h+,h-] , (2.18) 

we can write it as 

1 1 

S'efr[/l+, h'] = -hM^MNhN = 2^M(^A/Af + i^'^yMN)hN , (2.19) 



where 







&h-Sh- 

is the free theory differential operator acting on /i, and 



V= 5h+Sh+ 5h+&h- S2[h+,h-] (2.21) 



(5^ S_ 

,5h-5h+ 5h-5h 



is the part arising from the interaction with the matter fields cj). (Note that 52 also includes 
any counterterms proportional to the terms in the Einstein-Hilbert action which can appear 
for massive fields or when using Pauli-Villars regularization.) The functional integration in 
equation (2.12) leads to 



Z[J+,J"]~ j 5[/i+(r) - /l"(r)]e5''A^^Afiv^iv+iJiv/iivp/j ^ g-^Jiv/GMivJiv ^ (2.22) 
where G is the inverse operator of A, i.e. 

j Aam{x, y)GMBiy, x') d^y = 5'^{x - x')5ab ■ (2.23) 

Functionally differentiating the generating functional twice with respect to the source 
J (remember that Jm = (J^, — J^)), we get 



{m\VhA{x)hB{x')\m) 



= \Gab{x,x). (2.24) 

j=o 



\5Ja{x) i5JB{x' 

Specializing to the four possible values of the index pair AB, we get 

G++(x,x') = — i(in|T/i.(x)/i(x')|in) = Gf(x,x') (the Feynman propagator) 

G (x,x') = — i(in|T/i(x)/i(x')|in) = G-d{x,x') (the Dyson propagator) 

G ^(x,x') = — i(in|^(x)/i(x')|in) = G^(x, x') (the positive Wightman function) 

G-\ (x, x') = — i(in|/i(x')/i(x)|in) = G^{x,x) (the negative Wightman function)^ . 

In general, however, we cannot solve equation (2.23) exactly, and must use perturbation 
theory. Adopting a condensed matrix notation, in which equation (2.23) reads AG = 1, and 
assuming that we can invert A^ exactly to get the propagator G^ 

= 1 , (2.25) 



^Note that while the name Wightman is attached to those functions regardless of the order of perturbation 
theory, one usually uses Feynman and Dyson propagator only for the tree level result, while in higher orders 
they are known as time-ordered and anti-time-ordered or r functions. 
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we get from equation (2.23) (with appropriate boundary conditions) the exact recursion 
relation for G, 

G = G^ - K^G^VG . (2.26) 
Since ti^ is smah, we can solve perturbatively this recursion for G by iterating and truncating 

G = G° - K^G^VG^ + O {k^) , (2.27) 

which fully written reads 

Gab{x,x') = G%{x,x')-^^ jj G%dx,y)VMNiy,y')G%Biy',x')d\ d^y+O {k^) . (2.28) 

This is shown in terms of Feynman-Stiickelberg diagrams in figure 2. 

Since V is the second variation (2.21) of the S2[h'^ , h^] part of the effective action (2.18), 
the integrand in (2.28) is given effectively by S2[h'^,h^] with the fields h replaced by the 
zero-order propagators G^. We will calculate exactly this double integral, adding, of course, 
the tensor structure which we have ignored until now. 




Figure 2. The Feynman-Stiickelberg diagrammatic expansion corresponding to our calculation. The 
wiggly lines are gravitons and the straight lines are scalars. Summing the contribution of the loop, 
the tadpole (which vanishes in dimensional regularization for the massless, conformally coupled scalar 
in spatially flat FLRW backgrounds) and the counterterms, we obtain the (finite) dashed loop corre- 
sponding to V (2.21). 

It is worth mentioning that in the stochastic gravity framework [46-50] , one now rewrites 
the imaginary part of S2[h~^,h~] using an auxiliary stochastic field ^ to arrive at a real 
stochastic effective action, from which manifestly real and causal equations of motion for 
the perturbation h can be derived. An alternative expression for Gab can then be given 
as a sum of two terms called "intrinsic" and "induced" fluctuations, where the "induced" 
fluctuations are a direct consequence of the stochastic source. It is then important to consider 
the contributions of both terms, since otherwise secular terms occur in the two-point function 
for the metric perturbations. An explanation of this fact is given in section 4.4. 

However, in order to implement correctly the prescription for the adiabatic "in" vacuum 
in the interacting theory, it is more convenient to work directly with equation (2.28), as we 
will do here. 



3 The effective action for metric perturbations around a spatially flat 
FLRW background 

We consider massless and conformally coupled free scalar fields interacting with the metric 
perturbations around a spatially flat Friedmann-Lemaitre-Robertson- Walker (FLRW) back- 
ground driven by a cosmological constant. As described in the previous section, since we are 
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interested in calculating the correlation functions of the metric perturbations to leading order 
in in addition to the Einstein-Hilbert action for gravity with a cosmological constant 

we only need to include counterterms quadratic in the curvature for renormalization, 

S[g, 0] = ^ J [r- 2Aj d"x + Nai J [r''""' Rabcd - R'^^Rab 

I' ^ ^ I' 

+ Na2 R^y^gd^x-^Y. [{^""<Pk) (Vm<Pk) + mR<Pl 

k=i-' 

Note that we are considering this action in n dimensions so that we can use dimensional reg- 
ularization, and that all parameters are bare parameters, except for ^(n) = (n — 2)/[4(n — 1)] 
which depends on n such as to make the scalar field action conformally invariant in all di- 
mensions. In the large limit and after the rescaling of Newton's constant k? = Nk^, 
graviton loops are suppressed by higher powers of with respect to matter loops, and one 
can neglect them at leading order. The physical metric cjab is conformal to an almost flat 
metric Qab 

gab = a^{'n)gab = a^{v){^ab + Khab) , (3.2) 

with 7] denoting the conformal time. Neglecting graviton loops then amounts to truncating 
the effective action Ses, as given by equation (2.15), to second order in the metric perturba- 
tion hab- 

The functional integral over the matter fields was calculated in [38, 51] for the case of 
a single scalar field (A^ = !)• The basic idea is to use the conformal transformation formu- 
las given in appendix B together with the conformal invariance of the classical action for 
the conformal fields to rewrite equation (3.1) in terms of the metric Qab for small pertur- 
bations around flat space. One can then perform the path integral for the scalar fields in 
flat space treating perturbatively the interaction with the metric perturbations. In doing so, 
UV divergences arise which are canceled by an appropriate choice of the bare parameter (in 
dimensional regularization) 

Nai = Naf^'ifi) + a [(n - 4)"^ + In/i + O (n - 4)] , (3.3) 

where we have introduced the dimensionless constant 

a = ^ . (3.4) 

28807r2 ^ ' 

In contrast, the second bare parameter, Na^ = (a — /3)/12, does not get renormalized, which 
is a peculiarity of the conformal case [38, 52]. Taking all this into account, the final result for 
S'cff obtained in [38], which is valid to quadratic order in the metric perturbations, is given 

by 

S.^\tj = glb\ - 5 J + 3a 11 C+_(x)C7-«^-"(y)i^(x - y) d'x d^y 

+ // C+_(x)C+'^''-"(y)i^+(x-y;M)d^xd4y (3.5) 



2 
3 

2« 



C-,^Sx)C-^'"^\y)K- {x - y;f,) d'xd'y 
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with 



+ aj {R^^''''Rabcd-R''''Rmn) In 0^=5 

f / T~\ n 1 T 7 IT) T 7 \ 2 



12 



-gd X 



(3.6) 



+ 2ay" G""a-2(V„a) (V„a)7=^d^a; 

+ aja-^ (2aV™V„a - (V"a) (V^a)) (V"a) (V„a) 



"5 d^3; . 



where the divergences appearing in U and those of the bare parameters have canceled out and 
the remaining parameters are the renormahzed ones. Moreover, we have used the identity 



-K rtabcd 



nmn r> 
It rLr, 



-gd X 



i^abcd/^ 

^ ^abcdy 



-gd X, 



(3.7) 



which is vahd in four dimensions and follows from the Gaufi-Bonnet theorem and the defini- 
tion of the Weyl tensor in terms of the Riemann and Ricci tensors (A. 4). The kernels K are 
given by their Fourier transforms as 



K(x) 



K^{x;fi) 



-ITT 



&{-p')0{-p'')e 



0\Jpx 



dV 

("2^ 



In 



± i7r0(-p^ 



d^p 



(3. 



where ^ is the renormalization scale. It should be stressed that the dependence on ^ of these 
kernels and of the renormahzed parameters cancels out, so that S^s is invariant under the 
renormalization group. More specifically, we have 



A^a'r(/^o) -aln 



(3.9) 



for a fixed scale (Iq, and one can see that the logarithmic dependence on /x of the local terms 
quadratic in the Weyl tensor exactly matches that of the kernels K^{x; fi). Given any value 
of ai^^{^Q), one can always choose a renormalization scale /2 such that ai'^{fl) = 0. This 
is the choice that was made in [38], and we will also employ it throughout our calculation 
to deal with slightly more compact expressions. However, once we obtain the final result 
in the next section, we will specify the small changes needed so that the result is valid for 
arbitrary ^u. On the other hand, as mentioned above, the parameter 02 does not get divergent 
contributions for conformal fields and, thus, no dependence on /x. In this case the functional 
integration of the scalar fields just gives a finite contribution, —a/ 12, to the coefficient of 
the term with the square of the Ricci scalar. To take that into account, the new parameter 
P = a — 12Na2 has been introduced in equation (3.5). In [38] the particular choice (3 = a, 
corresponding to 02 = 0, was made. 

This effective action Ses in equation (3.5) should be understood as truncated to second 
order in hab, the perturbations of the metric gab around fiat space. The expansion of the 
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effective action up to that order can be obtained using tlie expansions around fiat space given 
in appendix A. 

In particular, from equations (A. 2) - (A. 4) we have 

Cabcd — ~^Tabcdmnpqd hF'^ + O (k ) 
'-^abcdmnpq Qr]p[aVb][cVd][nVq]m " QVm[aVb][cVd][nVq]p " QVm[ar]b]pVn[cVd]q + '^Tla[cVd\bTlm[pVn\q ■ 

(3.10) 

Taking into account that the scale factor a only depends on the conformal time r] (and 
denoting the derivative with respect to it by a prime), we then get 

+ iK^r.b,,™^"^'*^'^' j I [dmdnh^^{x)] [drdsK^ix')] K{x - x') d^X d^x' 

+ ^K^r,,,,™^'^'^*^'^' j j [dmdnh+,{x)\ [drdshl,{x')] K+ {x - x'; m) d^X d^x' 

- ^K^T.b.d^T''^'^'^^'^' j j [dmdnh„{x)\ [drdshl^ix')] {x - x';/x) d^xd^x' 
ST[9ab] = 4°^[a] + ^ab] + S^;?[a, Kb] , 

(3.11) 

where 

4°) [a] = ^ J [6aa" - 2Aa^ + at^a-^ {2a{a'fa" - {a'f) - 'ipK^a~\a" f] d^x 
s'q [a, hab] = - 

- ^/3ac j (2a2a(4) - Waa'a'" - 5a{a"f + 16(a')'a") h d^x (3.12) 
+ [aa" - 2ia'f - aK^a-^ {a{a'fa" - 2(a')')] d^x 
-pK f (^a^a^^) - Saa'a'" - a{a" f + U{a'fa"^ /i^o d^x 



2aa" - {a'f - Aa" + \aK^a-^ {b{a'f - Aa{a'fa") 



hd^x 



and Sq [a,hab\ contains all the terms which are quadratic in the metric perturbation hab- 
3.1 Semiclassical de Sitter background 

As explained in section 2.2, we will consider a background metric which satisfies the 
semiclassical Einstein equation (2.17). Assuming a FLRW background of the form a?{ri)rjab, 
the dynamics for the scale factor a{r]) is entirely determined by the Priedmann equation, 
which corresponds to the 00 component of the semiclassical equation. It can be immedi- 
ately obtained by functionally differentiating with respect to h?^ the expression for m 
equation (3.12). Taking into account that 5h{x) / 5h?^ {y) = —5'^^\x — y), we get 

6{a'f - 2Aa^ = 3aK^a-^{aY + 3pK^a~^ [2aa'a"' - a(a")^ - 4.{a'fa"] . (3.13) 

The quantum corrections on the right-hand side of equation (3.13) involve higher deriva- 
tive terms, which lead to spurious solutions that lie beyond the domain of validity of the 
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EFT expansion [45, 53]. A suitable way of dealing with that problem is the order re- 
duction method [54]. The basic idea, applied to this case, is to start with the equation 
6(a')^ — 2Aa^ = O (k^) and write a' in terms of a up to terms of order k"^, which is 



a' 



' = {2A/3)a^ + 0{k^). Differentiating this expression and proceeding recursively, one 
can express higher-order derivatives such as a" and a'" in terms of a up order k^. Next, one 
substitutes these results into the right-hand side of equation (3.13) and obtains an equation 
without higher-order derivative terms, but equivalent to equation (3.13) through order k^, 
which is anyway the order up to which that equation is valid within the EFT approach to 
perturbative quantum gravity. When doing so, the terms proportional to f3 cancel out. This 
was to be expected because they correspond to the trace of the functional derivative with 
respect to the metric of the term quadratic in the Ricci scalar, which is proportional to Di?, 
whereas the equation that we are substituting implies that R = A A, which is a constant. 
Taking all that into account, after employing the order reduction method we get 

3(a')^ -a(^ + ^an^A^ = O {k^) . (3.14) 

We can, therefore, define an effective cosmological constant 

yleff = a{i + ^aK^A)j , (3.15) 

which includes a small constant shift due to the quantum back-reaction of the matter fields 
on the background geometry. The scale factor is then given by 

a(r?) = --i- (3.16) 

with the Hubble parameter H = y^yleff/3. This corresponds to de Sitter spacetime, more 
precisely to the part which can be covered with a foliation involving flat spatial sections (the 
Poincare patch). 

We hereby rederive the well-known fact (see for instance [55-57]) that one can have 
self-consistent solutions of equation (2.17) which correspond to de Sitter spacetime with a 
modified cosmological constant, since, due to de Sitter invariance, the renormalized expecta- 
tion value of the stress tensor for the Bunch-Davies vacuum is a cosmological-constant-like 
term proportional to the metric. 

3.2 Tensor perturbations 

Concerning the metric perturbations around the mean background, we will from now on 
restrict to tensor perturbations (scalar and vector perturbations will be discussed in [58]), 
which fulfill the conditions 

/i = 0, /i°" = 0, 9^/i'"" = 0. (3.17) 

Inserting the background scale factor solution (3.16) into the gravitational action (3.12), 
the linear term Sq^ vanishes as expected. Integrating by parts and using 



/ 



h^'Vm'mnd'x = -\ I h^'^f'mrnnd^X, (3.18) 
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the quadratic term of the gravitational action becomes 



45jj^^[a, hab] = - (1 + (5a - 2P)k^H^) J a^{d'h"'''){dshmn) d^a; - QaK^H^ j a^h"''''h'^^ d^x 
+ 30^2 J {ah"'''){ahmn) In a d^x . 



(3.19) 



The effective action S^s (3.11) can then be written as 

Sesitb] = 4'^ Kb] - 4'^ Kb] + I [ Ph;,ix)) {n'h-P^{x')) K{x - x') d'x d'x' 



+ -^aK^ II (□/i+(x)) {n'h+P'^{x'))K+{x-x';fi)d^xd 



3 2 
-OK 

4 



//{□V,(.))C''.-"(.'))K-(.-.';P)dS.dV^ 



(3.20) 



Note that we can alternatively shift the last term in the gravitational action (3.19) into the 
kernels K^, which amounts to replacing fl with fia = —Jl/{Hri). Equation (3.20) corre- 
sponds to equation (2.18) and is the starting point for the next step: computing the metric 
correlations. 

It should be pointed out that when integrating out the matter fields, one naturally gets 
an effective action where, in contrast with equations (3.5) and (3.20), the differential operators 
are acting on the kernels rather than the metric perturbations [47] . If one assumes that the 
metric perturbations in the effective action are regular test functions, one can integrate 
by parts and recover those equations. This assumption is actually too restrictive when 
performing a calculation like ours and one gets boundary term contributions at the initial 
times, and t^, from the integration by parts of the time variables (the boundary terms 
at the final time T cancel out provided that it is larger than any other relevant times). 
Nevertheless, when applying the prescription described in section 2.1 and taking the limit 
— >■ — cx)(l ^ ie), such boundary terms vanish and do not contribute to our results. 



4 The gravitational two-point function 
4.1 Tree level 

To calculate the two-point function to first order in V (which is proportional to k?) as given 
by equation (2.28), we need to compute the lowest order propagator. For this we note that 
the gravitational action (3.19) is given at lowest order by 

= 0) = So = a\d'h^'^){dsh^n)d''x. (4.1) 

We now follow Ford and Parker [59] to express this action as that of two scalar fields cor- 
responding to the two physical polarizations of the tensor perturbations. We first Fourier 
transform the field with respect to the spatial coordinates, 

hmn{x) = j h^^{r^,pyP-^ (4.2) 
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and note that the gauge-fixing conditions (3.17) for tensor perturbations correspond to 

h = 0, /iO" = 0, p^/i"^" = 0. (4.3) 

Next, we choose a right-handed orthogonal set of three (real) vectors e^(p) such that 
p = Iple^ip), and define the following two (real) scalars corresponding to the two polar- 
izations: 

h+{v,p) = \h^n{ri.p) [eT{p)e1{p) - e^(p)e^(p)] 
hx{r],p) = hmniv,P)eTip)e2ip) ■ 

As is well known, such a definition is not possible in a way that would be continous for 
every orientation of p. However, with respect to the transformation p — )• —p one can adopt 
consistently the definition 

esi-p) = -esip) , ei(-p) = 62(2?) , e2{-p) = ei{p) . (4.5) 

Different definitions would be possible, but lead to the same result for the graviton propaga- 
tor (4.18). The relation (4.4) can now be inverted to give 

/i'""(r?,p) = -ih+{v,p) [eTip)enP) " e?(p)e?(p)] + h^iv^P) [eT{p)e^{p) + e^(p)e?(p)] • 

(4-6) 

Using the reality of hmn, i-e. hmniv^ ~P) = f^mniv^P)^ it follows that 



^0= ^ // a^iv) K{rj,p)h*;{ri,p)-p'h+{v,p)Kiv,P) 



di], 



(4.7) 



2tt 



+ h'^{ri,p)h*^{T],p)-p hyc{r],p)hl{r],p) 

which after inverting the Fourier transform reads 

So = -\j a" Wh+) (dM) + (d'h^) {d,hx)]d^x, (4.8) 

which is the action for a pair of minimally coupled scalar fields in a spatially flat FLRW 
background, whose canonically conjugate momenta are given by 

7r_|. = a^/i'_|_ , VTx = a^h'y. . (4.9) 

We quantize these fields by imposing that the field operators satisfy the canonical commu- 
tation relations at equal times, 

[h+{ri, x),TT+{r], y)] = [hy, (??, a;), ttx (r?, y)] = i5^{x - y) , (4.10) 

with all other commutators vanishing. Fourier transforming the spatial coordinates, this is 
equivalent to 

[/i+(7?,p),^+(??,q)] = [/ix(r/,p),^x(??,q)] =i(2^)^<5^(p + q). (4.11) 

From the action (4.8) it follows that the equations of motion for the operators in the Heisen- 
berg picture are 

- 2Hadr,) /!+=(□- 2Hadr,) hy, = . (4.12) 
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To calculate the propagator, we need to find a complete set of solutions ("mode func- 
tions") to those equations, which then determine the expansion of the perturbation operators 
/i+ and /ix and the associated vacuum state. In the following, we will perform the analysis for 
and trivially extend the results for hx- Fourier transforming the first equation in (4.12) 
with respect to spatial coordinates, it becomes 

~h'i{r],p) - ^h'4v,p)+p^~h+{7],p) = 0, (4.13) 

whose general solution is given by 

h+iv,p) = C{p)a+{p)i\p\rj - i)e-W^ + C*(-p)aV(-p)(|p|r? + i)e'IPl^ (4.14) 

with an arbitrary operator a+ and a normalization constant C, and where we have used the 
reality condition for h^{x), which corresponds to hermiticity for h^(x). For the canonical 
momentum, we get 

hiv,P) = --^^ [c(p)a+(p)e-'H'' - C*{-p)aU-py\P\^] • (4.15) 

The canonical commutation relations (4.11), together with the vanishing of the equal- 
time commutators for h-^- and tt-j- with themselves, are valid if and only if 

C{p)C*{p) = Ci-p)C*i-p) 

[a+{p),a+{q)] = [al{p) , al{q)] =0 

2\p\'C{p)C*{p) 

_ 1 

By choosing C{p) = H ^2|p|'^^ ^ , the operators a+ and fulfill the standard commutation 
relations for creation and annihilation operators. We then define the free field vacuum |0) by 

a+(p)|0) = ax(p)|0) = 0. (4.17) 

This vacuum which (for \p\r] ^ 1) resembles the Minkowski vacuum with respect to the 
choice of mode functions (4.14), is the Bunch-Davies vacuum [60] for the gravitons, and we 
will use it in the following (there is a subtelty in this case concerning IR divergences which 
is briefly discussed below). 

The lowest order two-point function of the tensor perturbations {Q\h°'^ {x^lf^ {y)\Q) can 
now be evaluated directly using the decomposition (4.6) into the + and x polarizations and 
the mode expansion (4.14). A straightforward calculation then leads to 



where 



and 



{G\h'^\n,pW\rj\qm = {27:f5\p + q)f{r^,r^\\p\)P-'-\p), (4.18) 

/(r/, r?', IpI) = -^(|p|r? - i){\p\rj' + i)e-'H(^-''') = /* (V, r/, |p|) , (4.19) 
2\p\ 



P^'^\p) = eUp)ei{p) + ei(p)e?(p) (el{p)ei{p) + e^,{p)ef{p) 



+ (eUp)e\{p) - e^(p)e^(p)) (eUp)ef{p) - e^(p)e^(p)) (4-20) 



padpbc _j_ pacpbd pabpcd 



- 18 - 



is the polarization tensor which can be written in terms of the projection tensor P"^ defined 

by 



•ab 



ab I s-ajrb 
1] i-0o"0 



p2 



(4.21) 



where = (0,p) in the coordinate system that we have chosen. The projection tensor P'^^ 
satisfies 



P^'Vab 



pab pcd^ 

P P ribc 



P' 



ad 



(4.22) 

From equation (4.18) we can now calculate all the lowest order CTP propagators for 
the graviton. They are given by 



GTdiV,V',P) = &''''M,p)Pabcd{p) , 

where the components of the 2x2 matrix are 

GO++(77,r?',p) = -i [0(77 - r]')f{r^,7^', \p\) + 0{rj' - v)nr^,r]', \p\)] 

G'+~{v,v',p) = -^r{v,v',\p\) 

G^~+{v,v',P) = -^f{v,v',\p\) 

(77, v',P) = -i [0{V - r]')nv, v' , \p\) + 0{r]' - r,)f{r], rj' , \p\)] 



(4.23) 



(4.24) 



Since / goes like |pp for small this two-point function must be truncated at some 
small momentum |p| = to have a well-defined inverse Fourier transform. The two-point 
function in coordinate space diverges then as In^ for ^ — t- 0, as was noted by a number 
of authors [61, 62]. This was to be expected since the same problem occurs for a massless 
minimally coupled scalar field in de Sitter space [63-65]. However, similarly to the stress 
tensor correlation function for a scalar field [65], gauge-invariant correlation functions of 
observables characterizing the local geometry in a finite-size region, such as the curvature 
tensor correlation functions, are finite and do not grow at large distances [65-67]. We will 
not discuss the issue in this paper, but plan to address it in a future publication [58]. 

4.2 One loop 

We can now use equation (2.28) to calculate the one-loop order graviton propagators which 
are of order in our perturbative scheme. As explained at the end of section 2.2, this is 
easily accomplished by replacing the fields hab in the effective action (3.20) with the lowest 
order propagators. Fourier transforming the spatial coordinates, we get 



GtUr].ri,V) = &TMri,v)-rPabcd{pWG'^^{r^,ri,p) 



\ABi 



(4.25) 



where 



G^^''(^,r/',p) = (5a-2/3)p2 / r-2G%,(77, r,p)eAfjvG^B(r, r/', p) dr 



+ (a + 2/3) y T-''{drG\ukn.r.P)yMN{drG\j,{T,rl,p))^T 



+ 3a 



(dl^v')G\u{^,T,v) 



{dl+p')G%B{r,r^\p) H-HT)dr 



-a 



(92+p2)G%,(r/,r,p)] VM^(r,r',p) [{dl + p')G%s{r' ,P)\ drdr', 

(4.26) 
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1 



^0 -1 
VMNir, t',p) 



and 



L{t -T',p;fi) P{t-t',p) 
-D{t-t',p) -L{t -T',p;fl) 



+ iiV(r-T',p) 



1 -1 
-1 1 



(4.27) 



This is written in terms of the Fourier transforms of the following real- valued kernels, in 
terms of which the complex-valued kernels in (3.8) can be decomposed: 



L{x;ii) 



In 



N{x) =71 0(-/)eiP^ 



d^p 
(2^ 



L{-x;n) 

: N{-X) 



(4.28) 



D{x) = ivrl 0(V)sgnpOei^'-^ = -D{-x) 



The Fourier transforms of those kernels can be calculated with the help of appendix D and 
are given by 



L{rj - r]\p) = cos [\p\{j] - r}')] V' 
D{rj — r]',p) = cos [\p\{r] — r/')] V 



1 

\ri — rj'l 
1 



7] — rj' 



(4.29) 



N{r] - r]',p) 



sm[\p\{r]-r]')] 
r] — r]' 



+ 7r6{r] - r]') 



After calculating the integrals using the ie prescription as explained in some detail in 
appendix E, we get for the two-point function of the metric perturbations 



{in\h-\r^,p)h^^{r^\q)\m) = {2TTf6\p + q)P 



jabmn , 



f{r],rj',p) [l + 6aK^H^ln 



{5a - 2/3) k'^H^ 



'^aK'^H'^(^h{'n,r]',p) - l2{ri,'n',p) - I^{r]' ,r],p) 
^aK^H^(l3{ri,ri',p) - h{r],r]',p) + h{r],r]' ,p) 



(4.30) 

where 

h{r],ri',p) = 2|prS'?'e"''^'^''~'''^ 

hir], r]',p) = |p|-3eiW(^+''') (|p|r, + i) (|p|r?' + i) [Fin (-2i|p|7?) + In (2i|p|?7) + 7] 

+ Ipp^'e-'H {\p\r] - i) (|p|?7' + i) In {-2\p\r]) 
hiri, V',P) = Ipr^'e-^IPK"-^') (IpI?? - i) {\pW + i) [in (2i|p|(7? - 7?')) + 7] 

hir], ri',p) = IpP^^'IpK^-'^') (|p|r? + i) (|p|r?' - i) [Fin (-2i|p|(7? - rj')) + In (2i|p|(r? - 7?')) + t] 
hiVir] ,p) = rj^{r]'f N{ri - r]',p) - iD{r] - r]',p) . 

(4.31) 
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This is the exact result for the one-loop-corrected (due to conformal matter fields) two-point 
function for tensor metric perturbations, and unlike previous results, it has been obtained 
for arbitrary pairs of times rather than just equal times. It is finite for non-coincident 
points and (apart from a global factor |pp^) invariant under the simultaneous rescaling 
{p — )• X~^p; 7], r\ — )• Ary, A??'}, which is a necessary (but not sufficient) condition for de Sitter 
invar iance. 

The origin of the two kernels N and D in the result for 15(7/, 7/',p) is the following. The 
two CTP propagators Gab in equation (4.26), or equivalently equation (E.3), appear with a 
linear operator acting on them which contains second-order time derivatives. If Gab is a true 

propagator rather than a Wightman function, i.e., for G++ and G (or G^^t in appendix E), 

the action of the linear operator gives rise to a term proportional to a Dirac b distribution in 
time. When that is the case for both G ab and Gcd in equations (4.26) or (E.3), the result of 
the two time integrals is simply a term proportional to the corresponding kernel, evaluated 
at the two times r/ and Vj in the argument of the metric correlation function. 

Equation (4.30) is valid for the particular choice of the renormalization scale such that 
a!^Y^{ji) = 0, as explained in section 3. For arbitrary values of the renormalization scale one 
should make the following replacement in equation (4.30): 



iVar(/i) + aln(|:) . (4.32) 



Note that this guarantees the invariance of the two-point function under the renormalization 
group, as follows from equation (3.9). 

4.3 The power spectrum 

In cosmology, two-point correlations are often expressed in terms of the power spectrum. To 
calculate it, we need to take the equal time limit r/' — )• r/ of the expectation value (4.30). 
However, this expectation value is not a regular function but a distribution with singular 
support, for which the equal-time limit does not make sense. Therefore, the concept of a 
power spectrum is, strictly speaking, ill-defined at the one-loop level (see appendix F for a 
more detailed discussion of this point). In order to compare with previous results in the 
literature, however, we can take the equal time limit of the regular part, which amounts to 
replacing the kernels appearing in I5 by 

AT/ I „^ ^ sin[|p|(r/-r/')] 
N{ri,r] ,p) - 



V-r]' (4.33) 
D{7]-T]',p) ^0. 

Following section 18 in reference [68] (and taking into account that our metric perturbations 
and theirs are related by a factor k), the power spectrum is defined to be a suitable normal- 
ization factor times the spatial Fourier transform of the contracted correlation function at 
equal times: 



S\\P\,V) = ^^IpfvamVbnj (in| /i'^^(r/, ^) /i™" (,?, 0) | in)e-'P- d^x 



(4.34) 
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which gives 



327r3 



2(1+ p^r?^) (^1 + 6aK^H^ In (^j^ 



(5a - 2/3) k'^H'^^ 



e2i|ph + ^g.^ (-2i|p|r/) + In (2i|p|7/) + 7] 



Uok'^H^ (1 + p%2) In (-2|p|?7) + Qan^H'^p^ri^ (2 - p^r?^) 



(4.35) 



Once again, for an arbitrary value of the renormahzation scale /i one needs the replace- 
ment (4.32). 

Let us consider two important limits of equation (4.35). It is useful to write it in terms 
of the physical momentum (of the physical spacetime gab)^ which is 



V 



a 



-H\p\r] . 



(4.36) 



The power spectrum depends on |p| and r/ only through the combination p, so that it is 
time-independent in the physical spacetime. 

For "sub-horizon" modes with p ^ H we should then recover flat space behavior, as 
the curvature of spacetime should have no effect on those modes. Taking therefore the limit 
H ^ 0, we get from equation (4.35) 



167r3 



p2 (1 _ 3aKV + O {k^)) + O (^j^ , (4.37) 



where we have used the asymptotic expansion of the Ein function from appendix C. We 
therefore have a power correction to the Minkowski spectrum [69]. 

On the other hand, for "super-horizon" modes with p <^ H we take the limit p — )• to 

find 



167r3 



1 + Qan'^H'^ (in + 7 j - (5a - 2^)k'^H^ + O (k^) 



+ 



(4.38) 



which corresponds to a small constant shift of the scale-invariant power spectrum for tensor 
perturbations in de Sitter obtained at tree level (for arbitrary renormahzation scale /i, one 
needs again the replacement (4.32)). This is in contrast with the logarithmic dependence 
on the comoving momentum |p| previously found for one-loop corrections to tensor [30] and 
scalar [28, 29] perturbations, but is in line with the more recent results on loop corrections 
for scalar perturbations due to minimally coupled fields obtained in [31], where the origin of 
this discrepancy was identified. 

4.4 Cancellation of secular terms 

In our result for the one-loop two-point function of the tensor metric perturbations, there is no 
dependence left on the initial times = to(l+ie) after taking the limit to —00. Being able 
to eliminate any dependence on to is crucial for obtaining a result which is invariant (aside 
from a global factor |p| '^) under the simultaneous rescaling {p — t- X^^p; r],r]' — t- Xr],Xr]'}, 
which is a necessary condition for the de Sitter invariance of our result. In contrast, a related 
study of one-loop corrections to tensor perturbations around de Sitter [27] found a result 
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which is particularly sensitive to large values of to and diverges in the limit to ~^ —oo. We 
provide here a simple explanation of the origin of those terms and why they are absent in 
our case. 

It is useful to consider the following decomposition for the anticommutator of the metric 
perturbations when the effects of matter loops are included: 

(4.39) 

where the noise kernel Nabcdix,y) corresponds to the symmetrized connected correlation 
function of the stress tensor 

Nabcd{x,x') = ^{0\{tk(.x),ta{x')}\0) - {0\fab{x)\0){0\ta{x')\0) . (4.40) 

In addition, h^i,{x) correspond to solutions of the linearized semiclassical Einstein equations 
^ab [9 ~^ — S^^x^ife^b + ^] / around a background metric Qab, {■ ■ ■)q denotes an average 

\ / ren 

over their initial conditions according to the initial quantum state of the metric perturbations, 
and G^l*^^ is the retarded propagator associated with the linearized semiclassical equation. 
Equation (4.39) can be derived using the stochastic gravity formalism, but also within a 
purely quantum field theoretical calculation [39, 46, 47, 49] and assumes a factorized initial 
state with the metric perturbations and the matter fields uncorrelated.^ The first term on 
the right-hand side of equation (4.39) can be interpreted as corresponding to the fluctuations 
of the metric due to fluctuations of their quantum initial state (including the effect of matter 
loops on their evolution as given by the linearized semiclassical equation) and are often 
referred to as intrinsic fluctuations. The second term corresponds to the metric fluctuations 
induced by the fluctuations of the matter fields. 

It turns out that the contributions from intrinsic and induced fluctuations exhibit secular 
terms when considered separately. The basic idea can be illustrated with the simple model 
of two bilinearly coupled harmonic oscillators [70, 71] with action 

S[Q, Q] = \j [mQ^ - MQ'^Q'^ + mq^ - muj'^q^ - 2cQq^ dry , (4.41) 

where tj denotes the time and overdots denote derivatives with respect to it. Here the 
oscillator q is the analogue of the matter fields, while the oscillator Q is the analogue of 
the metric perturbations. The analogous quantity to consider is the correlation function 
5^{0(^)) Q('?')}^ calculated perturbatively up to quadratic order in c. The free retarded 
propagator for the oscillator Q and the noise kernel characterizing the fluctuations of the 
oscillator q are given by 

MU ^442) 
N{r] - rj') = ^— cos [uj{r] - r/')] , 

®Such kinds of initial states give rise to certain pathologies [39, 70]. Hence, in our study we have employed 
the method described in section 2.1 to generate a properly dressed initial state for the interacting theory (with 
correlations between the metric perturbations and the matter fields). 
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where the noise kernel corresponds to the real part (the anticommutator) of the two-point 
function of the oscillator q in the interaction picture and we have assumed that it is initially 
in its ground state, so that ({g(^), 9(V)}) = (2ma;)~^ exp[ii7(r/ — r]')]. In that case the term 
for the induced fluctuations is proportional to 



The key point is that the sines and cosines involve exponentials with two opposite signs in 
the exponent. For resonant oscillators (with u = il) this will give rise to terms for which 
the T- or r'-dependent exponents (or both) cancel out and one is left with terms that do 
not oscillate as a function of r or r'. Therefore, when integrating over these variables, 
one obtains secular terms proportional to r] — to, r]' — Iq or both. The structure of the 
induced fluctuations for the loop diagram associated with the cubic gravitational interaction 
between the metric perturbation and the matter field is completely analogous, with a sum 
over the spatial momentum running in the loop, and for massless fields there is always some 
momentum for which the resonance condition is fulfilled. 

In fact, given a spatial Fourier mode p for the metric perturbations, the expression for 
the one-loop correction to its two-point function is the same as the contribution quadratic 
in the coupling constant c for a model consisting of an oscillator with i7 = |p| interacting 
independently (and bilinearly) with a set of oscillators {qj} with frequencies {u}j}. The result 
is a simple generalization of the two-oscillator model discussed above with an extra sum over 
the frequencies Uj. More precisely, one needs to integrate over a continuous w distribution, as 
dictated by the spatial Fourier transform of equations (4.28), which involves just an integral 
over and whose result is given by equations (4.29). In addition, to mimic the gravitational 
case for metric perturbations around de Sitter, one needs to consider the coupling between 
the oscillators to be time-dependent with terms proportional to ij. (One should also consider 
the frequency of the oscillator Q to be time dependent, but there the time-dependent terms 
are proportional to ry~^. Therefore, they are essentially irrelevant for the discussion about 
secular terms because they vanish as r/ — )• — oo.) 

Despite the existence of secular terms for the intrinsic and the induced fluctuations, they 
are absent in the result for the full correlation function because the secular terms from both 
contributions cancel out. Such cancellation can be inferred from our calculations and the 
final result in section 4.2, but can also be easily understood when considering the Feynman- 
Stiickelberg diagrams in the OTP formalism following the approach described in section 2.1. 
The basic idea is that for fixed arguments of the correlation function (associated with the 
external legs) and integrating over the times of internal vertices, all the propagators leaving 
a given vertex on the branch are of positive frequency^ with respect to the time variable 
of that vertex at sufficiently early times (the analogous statement with negative frequencies 
holds for vertices on the "— " branch). Therefore, instead of an equation like (4.43), one gets 
a sum of terms of the following form (or analogous forms): 



'^In the gravitational case this is a consequence of G"*" (rji , rj) being the Wightman function of the free theory 
associated with the Bunch-Davies vacuum, which involves an exponential factor with exponent — ii7(r;i — rj). 
In that case there are also factors involving powers of the conformal times, but the key role is still played by 
the oscillatory factors as far as the possible existence of secular terms is concerned. 




(4.43) 




(4.44) 
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where in this particular case all the r-dependent exponents have the same sign. This guar- 
antees that there will be no secular terms and the limit Iq — >• — oo will converge. It should be 
emphasized that within this formulation the terms that naturally appear, such as (4.44), al- 
ways contain a mixture of intrinsic and induced contributions. Note also that it is important 
to distinguish the initial times and for the two branches. Otherwise we would miss 
the exponential factor in terms proportional to e'*^'^"'"^''^*^~*o ), or e'*^'^"'"^''^*^~*o ) (appearing 
when one has a time-dependent coupling, like the cubic interaction terms in equation (3.1) — 
quadratic in the matter fields and linear in the metric perturbations — for a nontrivial FLRW 
background). Missing these exponential factors gives rise to spurious finite contributions or 
even divergences in the limit to ~^ "Oo for the second case. 

An alternative way of understanding the cancellation of the secular terms is by noticing 
that the calculation of correlation functions using the method described in section 2.1 gives 
the correlation function for the ground state of the interacting theory (or an adiabatic ground 
state for suitable time-dependent interactions). This is clearer for time-independent Hamil- 
tonians (including the interactions). The ground state is then a stationary state and the time 
dependence of the correlation functions only involves the relative differences between their 
arguments. Therefore, there cannot be secular terms (depending on some particular initial 
time) in the exact result, which implies a complete cancellation of secular terms order by 
order in the perturbative calculation. The lack of dependence on an arbitrary initial time and 
the corresponding need for a cancellation of the secular terms can be extended to the case of 
time-dependent interactions provided that appropriate adiabaticity conditions are satisfied. 

After the discussion in this subsection, the reason for the sensitivity to large values of 
to of the result in [27] is now clear: since only induced fluctuations were considered, there 
were secular terms analogous to those generated in equation (4.43) which did not cancel 
out. Note that in a similar calculation for flat space all the secular terms in the induced 
metric perturbations do cancel out. This is, however, an accident of such a special case:^ in 
general other cubic interaction terms (different from the gravitational interaction) will not 
lead to such a cancellation even in flat space. Neither is this a consequence of considering a 
sufficiently simple situation: there would be no cancellation for simpler non-derivative cubic 
couplings or even for the simplest case of two bilinearly coupled resonant harmonic oscillators 
with a time-independent coupling discussed above. 

5 Discussion 

In this paper we have obtained the exact result for the gravitational wave spectrum in the 
de Sitter spacetime including the one-loop corrections from a conformal scalar field. Our 
result is compatible with de Sitter invariance: besides spatial homogeneity and isotropy, it 
respects the isometry corresponding to a global rescaling of the spatial coordinates and the 
conformal time in spatially flat coordinates for the Poincare patch. In particular it contains 
no terms involving the logarithm of the co-moving momentum, which had been found in 
previous calculations of the one-loop correction to the power spectrum for both scalar [28, 29] 
and tensor perturbations [30]. Such explicit dependence on co-moving rather than physical 

®The reason for the cancellation when considering the gravitational case in a Minkowski background [48] is 
that the noise kernel for a massless field involves two global d'Alembertian operators which can be integrated 
by parts so that they act on the graviton propagator. The action of the Minkowski d'Alembertian on the 
graviton propagator gives a vanishing result in fiat space, but this does not happen in a general FLRW 
spacetime. 



- 25 - 



momentum is incompatible with the rescahng symmetry mentioned above, as pointed out 
in [31], and can lead to significant deviations from the standard picture for inflationary 
models lasting for a sufficiently long time. This problematic term is canceled out when the 
additional contribution from the spacetime volume measure is properly taken into account 
in dimensional regularization, which gives rise to a term proportional to the logarithm of the 
scale factor in equation (3.5). The cancellation was explicitly shown for scalar perturbations 
in [31], but although expected, it had not been shown for tensor perturbations until now. 
This is, in fact, a rather generic phenomenon for one- loop calculations involving light fields in 
de Sitter and a similar cancellation in a slightly different context had also been emphasized 
in [72]. 

In contrast with existing computations of one-loop corrections to cosmological pertur- 
bations, where only equal times were considered, we have obtained the correlation function 
for arbitrary values of the two times. This is important so that one can check the full de 
Sitter invariance of the correlation function (the condition mentioned in the previous para- 
graph is a necessary but not sufficient condition for de Sitter invariance) since general de 
Sitter transformations do not preserve the equal-time condition. Furthermore, whereas cal- 
culations so far have only been approximate and typically valid for sufficiently late times, so 
that the co-moving mode under consideration is well outside the horizon, our calculation is 
exact and valid for arbitrary times, which makes it possible to check exactly whether de Sitter 
invariance is preserved by one-loop corrections in our case. Moreover, since we considered 
arbitrary pairs of times and explicitly performed the regularization and renormalization of 
the UV divergences arising in our one-loop calculation by using dimensional regularization 
and introducing suitable local counter-terms in the bare gravitational action, we could ex- 
plicitly see that although the renormalized two-point function is a well defined distribution, 
the coincident-time limit in spatial Fourier space diverges. (In position space, however, the 
equal-time limit is finite for spatially separate points.) That means that, strictly speaking, 
when loop corrections are included, the power spectrum is no longer a well-defined quantity 
in general.^ A more detailed discussion on this subtle point is provided in appendix F. 

Some earlier calculations of loop corrections from conformal fields to cosmological per- 
turbations (see [27] and references therein) found terms in the result for the two-point function 
which are proportional to the initial time to and diverge as one takes the limit Iq — )• — oo. This 
would imply a bound on the number of e-foldings of inflation, established by requiring that 
the amplitude of primordial cosmological perturbations does not exceed the limits imposed 
by CMB and large-scale structure observations [27]. Moreover, the need to consider a finite 
value of to would prevent de Sitter invariance since the terms with explicit dependence on 
to are incompatible, for fixed to, with the rescaling symmetry mentioned above (in our cal- 
culation instead any terms with dependence on to disappear when properly taking the limit 
to — )• — oo). This question was discussed in detail in section 4.4, where these problematic 
terms were identified as secular terms that arise due to the cancellation of the exponents of 
oscillatory factors when calculating the partial contribution to the two-point function often 
referred to as induced fluctuations. Furthermore, we also explained on general grounds why 
such cancellation of oscillatory factors and the corresponding secular terms should not be 
present in the result of a perturbative CTP calculation of the full correlation function, a 
point confirmed at one-loop by our explicit result of section 4.2. Finally, we clarified that 
considering only the contribution of the induced fluctuations without also taking into account 

®This point had not been appreciated in previous calculations because all UV divergences in the power 
spectrum were removed by hand rather than through a detailed process of regularization and renormalization. 
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the one-loop correction to the evolution of the intrinsic fluctuations is the crucial reason for 
the lack of cancellation of the secular terms in [27]. It should be pointed out that other quan- 
tities (different from quantum correlation functions calculated using the approach described 
in section 2) can exhibit secular terms. However, that does not necessarily mean that those 
quantities can grow significantly: it may just signal a breakdown of perturbation theory for 
sufficiently long times, while a non-perturbative calculation would reveal a perfectly regular 
behavior. This situation is illustrated by the analysis in [73], where the evolution for ar- 
bitrary long times of linear perturbations of the semiclassical background around de Sitter 
(including effects of matter loops) is obtained by solving non-perturbatively the equation 
governing their dynamics. 

We close this section briefiy describing our plans for future work in order to extract fur- 
ther relevant information from the results obtained here and to generalize them. One of our 
goals is to calculate the one- loop Riemann correlation function [58], which includes contri- 
butions from scalar, vectorial and tensorial metric perturbations. The tensorial contribution 
is the technically and conceptually most elaborate one, and the results obtained here are the 
key ingredient for achieving our goal. The Riemann correlation function (with appropriately 
raised indices) exhibits a number of interesting features. First, it is gauge-invariant even 
when the corrections due to matter loops (but not to graviton loops) are included. Second, 
it is compatible with de Sitter invariance (in contrast with other gauge-invariant quantities 
more commonly employed in cosmology). In fact, whether de Sitter invariance is preserved 
can be manifestly seen by expressing the result in terms of maximally symmetric bitensors 
[74, 75]. So far the rescaling symmetry of our result mentioned above is a necessary but not 
sufficient condition for de Sitter invariance since it only corresponds (together with the E(3) 
symmetry of the spatially flat sections) to a subgroup of the full isometry group, 0(4,1). 
Moreover, the additional isometrics act nontrivially on the tensorial structure and the use 
of maximally symmetric bitensors becomes particularly convenient. Thirdly, the Riemann 
tensor provides a complete and suitable characterization of the local geometry (geometrical 
observables associated with a finite region of spacetime can be obtained from it) while being 
insensitive to potentially IR divergent super-horizon contributions which cannot be probed 
by observables involving finite spacetime regions. 

A particularly interesting generalization of our work is to consider massless minimally 
coupled fields instead of conformal ones since they typically give rise to much larger IR effects 
in de Sitter, and can even provide some hints on the viability of the proposed mechanism 
for a secular screening of the cosmological constant in the pure gravity case [2, 3]. Having 
considered the conformal case first has, nevertheless, provided a good testing ground for our 
techniques. Furthermore, since it is a UV issue, the conclusion that the one-loop correction 
to the power spectrum is strictly speaking divergent, as discussed in appendix F, applies 
in general (including the minimal coupling case). Similarly the cancellation of the secular 
terms discussed in section 4.2 also holds generically because its main contribution comes from 
modes well within the horizon. In this respect, an additional contribution of our work is to 
shed light on the claim of [27] that even loops of conformal fields can give rise to large effects 
on cosmological perturbations. 
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A Metric expansion 

For the perturbed metric gab, we have 
gab = Vab + i^Kb 

gab ^ ^ab _ ^^ab ^ ^^^a^^bm ^ ^^3^ 

^ = 1 + L/j + L^h" - -yhmnK^'' + O (/l^) . 
Z o 4 

where = IGttGn and all indices are raised and lowered with the unperturbed metric "qab- 
For the calculation of the Christoffel symbols (the connection) and the curvature tensors 
we will regard hab as a tensor field in fiat space. Moreover, we use □ = rf^^dmdn and 
Qm _ ^mnQ^ Taking that into account, for the Christoffel symbols we get 

n"c = l^s\, - In^K^s^-,, + o [h^] 

S ac — dch 

The calculation of the curvature tensors is then a simple exercise, 

Rab = \k {d,nS"\b - dadbh) - h^^dy,,S"\-^a " \^^'' VmnV^^ ^y.S^ b]a + O {h^) 

R = K id^dnh""'' - ah) + {dndmh + - 2dnd''hma) (A.3) 

Finally the Weyl tensor is given in four dimensions by 

Cabcd = Rabcd — ^ i^acdbd — Raddbc + Rbd9ac — Rbc9ad) + -^R{9ac9bd — 9ad9bc) = ^ (h) ■ (A. 4) 

B Conformal transformation 

Under the conformal transformation 

9ab = e^'^9ab, (B.l) 



(A.2) 
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Jb - gbcg""") dmU} 



the transformed Christoffel symbols are given by 

and the curvature tensors foUow as 

Rbd = Rbd-'2 [VbVdOO - {VbUj){VdUj) + 5M(V™L^)(V„tj)] - 9mV™V„w 
e^'^fi = R- eV^VmUJ - 6{V"'uj){VmUj) , 

where the covariant derivative V refers to the metric gab- 



(B.2) 



(B.3) 



OO JU 



C Special functions 

We define the entire function Ein(z) by 

Qt - I ^ 

Ein(z) = / dt = > , , , 

Its asymptotic expansion at infinity (r — t- oo) is given by 

Ein (ar + (3) 7 - hi (-(ar + /3)) + e""'+^ 



1 1-/3 



(C.l) 



(C.2) 



where 7 is the Euler-Mascheroni constant. 

The following useful integrals involving Ein(z) can be easily obtained by partial inte- 
gration: 



g"^* 1 ac 

dt = -e" b 

bt + c b 



Ein 



(bt + c) +ln(6t + c) 



e"* ln(6t + c) dt 



e«* -e~ b In {bt + c) - e~ b Ein - {bt + c) 



1 ac 

e"* Ein(6t + c) dt = -e"' Ein(6t + c) + -e" b 
a a 



Ein ( ^ (^i + c] 



„ /a + ft /, 
Ein I — — {bt + c] 



(C.3) 



These integrals depend continuously on the parameters a and c, so that for instance the 
integral of Ein(5t + c) can be calculated by taking the limit a — t- on the right side of the 
last integral. 



D Distributions and their Fourier transforms 

Here we analyze the main distributions which appear in section 4.2. Good references for the 
general theory of distributions are [76] and [77] . 

The principal value distribution, denoted by P - , is defined as 



t e^O t 



(D.l) 
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where d-lim stands for "limit in the sense of distributions" , i.e. 

e{t - e) + 0{-t - e) 



f{t)dt 



(D.2) 



where f{t) is an appropriate (fast decaying) test function. The Fourier transform can be 
calculated to give 



1 



,iLjt , 









lim 


[(/>r)"H 


= lim 






e->0 



lim / dt + lim / dt 

) — >^oo /_„ t r— >oo / t 



lim 



liin [Ein {iojt) + In + Um^ [Ein (iwi) + In 
= lim [Ein (— iwe) + ln(iL<j) — ln(— iw) — Ein (iwe 
iarg(ia;) — iarg(— icj) =i7rsgnci;. 



(D.3) 



Another distribution of interest is V' j which is defined as 



P«= d-lim 



€^0 



0{t - e) 



+ (5(i)(ln(/xoe)+7) 



(D.4) 



where /xq > is an arbitrary reference energy scale to make the definition dimensionally 
correct, and we use the symbol V to distinguish from the previous V which did not include 
a S distribution. Its Fourier transform can be calculated by proceeding as in the previous 
case, and is given by 



J r'^e''''dt = lim 



e->0 



dt + ln(//oe) + 7 



= lim lim [[—7 — ln(— iw) — In r] + In r — Ein (iwe) — In e + In e + 7 + In //q] 

e— >0 



= ln 



\lo\ 



ITT 



(D.5) 



Using those two results we may write the following compact expression: 
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d-lim 

e-)-0 



t t 

\t\ 



+ 2,5(t)(ln(Mo€)+7) 



(D.6) 



where we denote the distribution defined by the last line as T^'ltp^. Simply applying to this 
equation the formula for the (inverse) Fourier transform of a function with the argument 
shifted by a constant, one gets 
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We have now all the elements for calculating the temporal inverse Fourier transform of 
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First, from equation (D.7) we obtain 
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cos(|p|t)P'^-<5(t)lnf^ 



Second, writing 0[—p^) as 1 — (9(|p|^ — we can immediately calculate 

Combining the two we finally get 
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which provides the kernels (4.29) which we need for our calculation. 
E Calculation of the integrals in equation (4.26) 

Here we will explain in some detail how to calculate the integrals in equation (4.25) using 
the ie prescription described in section 2.1. 

The first thing that we will do is to substitute the definitions of the distributions (4.29) 
from appendix D into these integrals. We then note that the (5(t)ln(^o) part from equa- 
tion (D.4) appearing in the kernel L can be used to integrate over r', so that the logarithm 
can be combined with the one appearing in the third integral of equation (4.25). Setting 
/io = /U, this gives 



(9,^+p2)G0(r,,r,p 



1 
-1 



ln(|:)dr. 
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We will perform this integral explicitly, which will serve as an illustration of the computation 
of the other integrals since they all follow the same general pattern. 

The two-point function that we will calculate is the imaginary unit times the positive 
Wightman function G ^ defined in (2.24). The contribution of the above integral is then 
given by 

(E.2) 
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We are integrating from tlie initial time to the final time T, according to whether the 
integration variable r belongs to the upper or lower half of the contour shown in figure 1 . 
For the Wightman functions we have 



(E.4) 



since these functions are proportional to / or /*, which are solutions of the equation of 
motion (4.13). To make sure that we do not get a contribution from the arbitrary final time 
T, we first write 
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&_{r,,r^\p) = -GU7^,v',p) + G\irj,7^',p), 

which simplifies J to 
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We see that the result will not depend on T provided that T > rj, r]' , as expected. Next, we 
calculate 
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and together with equation (E.4) this leads to 
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These integrals can now be easily computed. Using the definition of the Wightman functions 
from equation (4.24), we obtain 
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Finally, we take the limit — )• — oo(l =F i^)) which causes the exponentials containing to 
vanish. After collecting terms we get 

J = i^ (Ipir, - i) (Iplry' + i) e^l^l = 2i/(r/, r,\p) , (E.IO) 

so that the contribution of the integral (E.l) is 

6aK2ln(^)/(r/,7?',p) (E.ll) 

as given in the final result (4.30). 

For the remaining integrals containing the kernels (4.28), we need the special function 
Ein, which is given (together with the necessary indefinite integrals) in appendix C. To further 
facilitate the calculation, it is convenient to combine the integrals containing the kernels L 
and as well as D and N to give the distributions 

Q^M{r-T')^_L_ ^ (E.12) 
r — T 

which allows to take the limit — )• — oo(l =p ie) already in the result of the integral over r' 
before integrating over r. 



F Coincidence-limit divergence 

The one-loop two-point function of the metric perturbations obtained in section 4.2 is a 
well-defined spacetime distribution, which can be obtained by computing the inverse Fourier 
transform of equation (4.25). This means that when integrated with appropriate spacetime 
test functions, it will give a finite result. However, if one takes the equal-time limit ry' — t- 
the resulting two-point function is no longer well defined as a distribution^*^ and its spatial 
Fourier transform, which corresponds to the power spectrum, is divergent. It should be 
emphasized that the usual divergences arising in our one-loop calculation have already been 

^^Mathematically, this corresponds to the fact that given a distribution in a manifold, its restriction to 
a submanifold is in general not guaranteed to be a well-defined distribution. The restriction should be 
implemented by considering a sequence of test functions on the full manifold with decreasing width around 
the submanifold, applying the distribution to the test functions and seeing if the limit of the resulting sequence 
(the zero-width limit) is well-defined. 
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properly regulated and canceled out by local counterterms in the bare gravitational action. 
There is nothing else to be subtracted: the two-point function is a well-defined distribution 
and it will give finite results when employed to calculate proper physical observables. It 
is just that when including loop corrections, the power spectrum is strictly speaking not a 
well-defined observable any more. 

We elaborate below on this point, which might sound surprising at first, and illustrate 
it with a simple example in flat space that shares all the key features. The fact that this 
issue also appears in flat space is not unexpected since it is a UV effect. 

F.l Simple example in fiat space 

Let us consider the following connected two-point function for a massless real held (j){x) in 
the Minkowski vacuum: 

2_N(x,x') = l{{^\x)J\x')})^, (F.l) 

where we introduced the notation (^A{x)B(x')'^ = (^A{x)B{x')'^ — (^A{x)'^(^B{x')'^ . It cor- 
responds to the real part of the product of two Wightman functions — 25Re x')]^. One 
can calculate it in (full) Fourier space as a convolution over four-momentum [48], and the 
result is 

N{p) = 7r0(-/) . (F.2) 

Note that we included a factor (Svr^)"^ in equation (F.l), so that the kernel N{x) coincides 
with the one introduced in equations (4.28). The expression in spacetime coordinates is then 
obtained by calculating the inverse Fourier transform and is given by 

N{x-x') = \ v ^, (F.3) 

vr2 ((^_;^')2)2' ^ 

where V denotes Hadamard's flnite part distribution [76]. Hence, we can explicitly see 
that N(x — x') is a well-deflned spacetime distribution and its Fourier transform is finite. 
Alternatively one can consider the spatial Fourier transform 

Ni, -,',p) = -^MpM^ + _ ,0 , (F.4) 
r] — rj' 

which is again a well-defined distribution, but diverges in the equal time limit rj' — )• ij, since 
it has singular support there. Hence, the power spectrum associated to N{x,x') is divergent. 
In contrast, the Fourier transform with respect to time is given by 

Nioj, x-x') = 8n^% ( - cos(L.r)^ , (F.5) 

where r = \x — x'\, and its spatial coincidence limit N{uj, 0) = is regular. 

There is, in fact, a relatively simple qualitative explanation for the divergence of the 
power spectrum associated with the product of two Wightman functions as well as for the 
radically different behavior of the temporal and spatial Fourier transforms [78]. The diver- 
gence of the equal time limit for the spatial Fourier transform (and thus the power spectrum) 
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can be easily understood if one writes the product of Wightman functions as a convolution 
in 3- momentum space: 

G+{v-v',P-q)G+{v-v',q)^,=- / 7^, (F.6) 

where = For rj = rj' power counting immediately reveals that the integral diverges. 
Note that momenta with arbitrarily large modulus \q\ contribute to the momentum integral 
even when keeping the total momentum p fixed, because they appear with opposite signs in 
the two Wightman functions. This is crucial to understanding the different behavior of the 
spatial and temporal Fourier transforms: although the Wightman function only involves pos- 
itive frequencies, it admits both positive and negative spatial momenta along any direction. 
More specifically, one can write the product of the two Wightman functions as a convolution 
in full Fourier space^^ 



1 



^ipO {ri-r]') ^-ip{x-x') 



(F.7) 



where q^ = u)q = \q\ in the last equality. The spatial coincidence limit of the temporal Fourier 
transform corresponds to fixing p^ and integrating over p the last integral in equation (F.7). 
This double integral over p and q is finite because \q\ is bounded by the Heaviside function, 
which corresponds to the positive frequency condition of the Wightman function and implies 
< In contrast, the equal time limit of the spatial Fourier transform amounts to fixing 
p and integrating over p^ the last integral in equation (F.7), which essentially eliminates 
the bound on |q| imposed by the Heaviside function and leads to a divergent result due to 
contributions with arbitrarily large \q\. In fact, the result of integrating over p^ corresponds 
exactly to the right-hand side of equation (F.6). 

We have considered the example of the correlation function (F.l) in fiat space for sim- 
plicity, but the essential points also hold for other correlation functions calculated at one-loop 
order in flat space. In fact, the anticommutator of the stress tensor can be obtained by apply- 
ing differential operators to (F.l) [48]. For a conformal field its projection to the TT sector 
is given by 



^({T.,(x),f,,(x')})^^ = ^/ e^^(— v.,,, {d', + \p\'yN{r,-r^',p)^, (F. 



where the TT projector Pabcd has been defined in equation (4.20). Equation (F.8) can be 
immediately extended to a n-dimensional spatially fiat FLRW spacetime with scale factor 
a(?7) if one multiplies the fiat space result by a?'~^{ri) o?~^{r]') [47, 79]. Moreover, since it 
is a UV effect, one expects that the main conclusions in this appendix will also apply to 
Hadamard states in general curved spacetimes. 

It should be stressed that equation (F.8) enables one to anticipate easily the singular 
coincidence limit of the power spectrum for the tensor metric perturbations even before cal- 
culating the time integrals corresponding to the interaction vertices in the Feynman diagram 



^^Equations (F.6) and (F.7) can be easily generalized to the massive case by taking cjp = ^/p^ + m? and 
adding rn^ to the arguments of the delta functions in equation (F.7), so that for instance becomes + m? . 
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for the one-loop correction of the metric two-point function, as carried out in section 4.2. 
The argument is the fohowing. First, one takes into account that the symmetrized connected 
two-point function of the Einstein tensor is given at one loop by the symmetrized connected 
two-point function of the stress tensor in the de Sitter background (the amputated loop 
diagram) [75]. Note that such a two-point function is given by a well-defined spacetime dis- 
tribution and it is not affected by the contributions from the renormalization counterterms. 
As mentioned in the previous paragraph, for conformal fields the stress tensor two-point 
function in a spatially flat FLRW metric is directly related to the result in Minkowski space- 
time, whose projection to the tensorial part is given by equation (F.8). One can see that it 
contains terms with up to four time derivatives acting on the Dirac 5 distribution function 
contributing to the kernel N(r] — i]',p), as given by equation (F.4). Since the two-point 
function of the Einstein tensor results from applying a fourth-order differential operator to 
the two-point function of the metric perturbations, one can conclude that the spatial Fourier 
transform of the latter should necessarily contain terms which diverge like a 6 distribution 
in the equal-time limit. 

One might try to circumvent the singular coincidence limit by defining the power spec- 
trum not as in equation (4.34), but as a constant times the quantum-corrected graviton mode 
functions derived from the semiclassical Einstein equations, along the lines of the treatment 
of Park and Woodard [80, 81] for the massless, minimally coupled scalar. 

F.2 Positivity of the power spectra 

Given a Hermitian operator Q, its two-point function must be positive semidefinite since 
QQ) = (Q^Q) > 0. The same is true for the connected two-point function (^Q^^ 

because for any given state Q — is also a Hermitian operator. In particular, one can 

consider the operator (/)^(x) smeared with a suitable (rapidly decreasing) test function F[x) 
as an example of such Hermitian operators: 

Q = f F{x)^'^{x)<l^x. (F.9) 



(Note that although (j)^{x) is a divergent composite operator, (f)'^{x) — (^(jP'{x)j is finite.) From 

the argument above one can immediately conclude that the kernel N{x, x') in equation (F.l) 
is positive semidefinite in the following sense: 



S[F] = JJ F{x)N{x,y)Fiy)d^xd'y = jj /,(p)|'/,o ^',P)/.o(^') dr dr' > , 

(F.IO) 

where we have assumed spatial translation invariance and that the test function factorizes 
into a temporal and a spatial part: F{x) = fs{x)frig{rj), with frjoit]) being some test function 
centered around r/o, such as 



1 



U{v) = ^^e . (F.ll) 

If we took the limit at — s- 0, we would be naively led to conclude that the associated power 
spectrum 6p{rj) ~ N [rj, r];p) is positive semidefinite. 

However, if one takes the simple flat space example given by equation (F.l) and neglects 
the local term, which diverges in the equal time limit, one gets Sp{rj) ~ N [rj, r],p) = —\p\ < 0. 
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This highlights the importance of the local term, which cannot be neglected in general. 
Indeed, if we keep the local term, we have 



S[F] 



lJUr} 



sin[|p|(T-r')] 



^(r')drdr' + 7r 
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(F.12) 



where we considered the limit \p\crt <C 1 in the last equality. 

Similar conclusions apply in general to other two-point functions at one loop. An 
immediate example is the connected two-point function of the stress tensor for the Minkowski 
vacuum, whose transverse traceless projection is given by equation (F.8) and to which the 
previous arguments can be straightforwardly extended by smearing with tensorial test fields 
F"'^{x) instead of test functions F{x). Furthermore, as already mentioned above, for a 
conformal field the connected stress tensor correlation function in n-dimensional flat space is 
simply related by a factor a?~'^{'q) a'^~'^{r]') to the result in an arbitrary spatially flat FLRW 
spacetime with scale factor a{rj). The correlation function of the metric perturbations, which 
we have calculated here perturbatively up to order for tensor perturbations around de 
Sitter would be another example. In this case, the positivity requirement only implies a 
positive result for the tree-level contribution (of order k^) in general, since the perturbative 
expansion is only asymptotic, i.e. valid for k — )• 0. Nevertheless, if one considers correlation 
functions for other quantities that can be obtained from the metric two-point function but 
for which the tree level contribution vanishes, positivity enforces a nontrivial condition for 
the loop correction (of order k^) as the first term in this asymptotic expansion. One such 
example is the connected two-point function of the Ricci tensor (with one index raised), 
which vanishes at tree level for perturbations around a de Sitter background, and essentially 
coincides with the connected stress tensor correlation function. 

F.3 Physical implications 

In section 4.2 we found a divergent result for the renormalized one-loop correction to the 
power spectrum of tensor metric perturbations. As illustrated with simpler examples and 
explained with general arguments in the previous two subsections, this is in fact a fairly 
generic situation for one-loop corrections. Here we elaborate on the physical interpretation 
and implications of this fact. 

As emphasized above, the spatial Fourier transform of the two-point function is a well- 
defined distribution in terms of the two times: it is only the coincidence limit which is 
divergent. Truly physical observables, nevertheless, will involve integrals over the two time 
arguments and will hence be finite. Furthermore, although the underlying reason for the 
divergence of the power spectrum is of UV nature, involving points with small invariant 
separations, it can easily "contaminate" attempts to extract IR features if one is not care- 
ful enough. In particular, simply considering spatial Fourier modes with sufficiently long 
wavelength does not necessarily exclude UV effects (unless one smears over sufficiently long 
times as well): calculating the spatial Fourier transform for equal times involves integrating 
over points with arbitrarily small spacetime separations. Alternatively, this can be seen in 
Fourier space as a consequence of modes with spatial momenta of arbitrarily large magni- 
tude contributing to the loop. Their contribution can only be suppressed by suppressing 
high frequencies, which can be achieved by considering sufficiently large time separations or 
integrating over sufficiently long times. 
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One way of extracting IR properties even in the equal-time limit would be to consider 
correlations between large spatial regions with no overlap. More specifically, suppose that 
we smear the two-point function with two different smearing functions Fi{x) and F2{x) 
with no overlapping spatial support. The corresponding smeared correlation function will 
be finite even when at — t- 0. The relation to actual measurements, however, deserves further 
discussion. Let us consider observables with a certain spatial extension and time duration 
at, characterized by smearing functions Fi{x) and F2{x). One can think of the ensemble 
with the outcomes of measuring these observables if they were repeated many times. For 
very small at the statistical dispersion for the measurements of the observables associated 
with Fi{x) or F2{x) would become very large even if their spatial support is large. On the 
other hand, the correlator of the two would be finite even in the limit at — > if the spatial 
supports of Fi{x) and ^2(2;) do not overlap: each one exhibits large fluctuations, but they 
are mutually uncorrelated. Note, however, that the product of the two measurements will in 
general be very large for each realization and only the average over many realizations would 
agree with the quantum correlation function and have a moderate value. 

Let us illustrate all this with a concrete example: the current measurement of CMB 
anisotropics. If one were to perform the measurements for sufficiently short times, loop cor- 
rections from QED would eventually give rise to large fluctuations, which would be due to 
UV effects of the current vacuum polarization, having nothing to do with IR effects of cos- 
mological origin. Such large statistical fluctuations would go away, after averaging over the 
outcomes of a large number of repeated measurements, if one correlated the measurements 
corresponding to non-overlapping solid angles with non-vanishing angular size. Neverthe- 
less, such spurious statistical fluctuations can only be avoided for single realizations if the 
measurement lasts for a minimum period of time. 
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